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The disciplines of asteroseismology and extrasolar planet science overlap methodically in the branch of high-precision 
photometric time series observations. Light curves are, amongst others, useful to measure intrinsic stellar variability due 
to oscillations, as well as to discover and characterize those extrasolar planets that transit in front of their host stars, peri- 
odically causing shallow dips in the observed brightness. Both fields ultimately derive fundamental parameters of stellar 
and planetary objects, allowing to study for example the physics of various classes of pulsating stars, or the variety of plan- 
etary systems, in the overall context of stellar and planetary system formation and evolution. Both methods typically also 
require extensive spectroscopic follow-up to fully explore the dynamic characteristics of the processes under investigation. 
In particularly interesting cases, a combination of observed pulsations and signatures of a planet allows to characterize 
a system's components to a very high degree of completeness by combining complementary information. The planning 
of the relevant space missions has consequently converged with respect to science cases, where at the outset there was 
primarily a coincidence in instrumentation and techniques. Whether space- or ground-based, a specific type of stellar pul- 
sations can themselves be used in an innovative way to search for extrasolar planets. Results from this additional method 
at the interface of stellar pulsation studies and exoplanet hunts in a beyond-mainstream area are presented. 
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1 Introduction 



The (transiting) extrasolar planet fields and the asteroseis- 
mology field see a convergence of instrumentation that cul- 
minates in the insight that beyond this purely technical level, 
a much more fundamental connection exists in the shared 
desire for the most exhaustive characterization of stellar and 
planetary systems at all possible with the available diagnos- 
tics. I will first sum up the relevant current context with an 
emphasis on planets around evolved stars, and then specif- 
ically address the topic of oscillation timing as a means to 
detect planets. 

In terms of successes to detect planets around evolved 
stars, the subdwarf B stars as host stars stand out as a group. 
This class of evolved objects will be introduced, and the dif- 
ficulties in explaining their sheer existence mentioned. Inci- 
dentally, the asteroseismology of subdwarf B stars is also a 
very active field. The EXOTIME planet searching program 
will be presented, which takes advantage of the long-term 
behaviour of these pulsations. Current ideas on subdwarf B 
evolution, and the potentially crucial role of planets, will be 
discussed, frequently resorting to the V391 Pegasi system. 



2 Extrasolar planet detection methods 



2.1 Overview 



The first extrasolar planet candidate was observed in ll989L 
without however b eing claimed as such at the time. 
Latham et alj (1989) instead suggested that the sub-stellar 
companion they had found around HD 114762 probably 
was a brown dwarf: a class of intensely searched-for objects, 
yet mostly elusive, at that period, and found to be intrinsi- 
cally rare in the role of companions to normal stars today. 

The first detection of an extrasolar planet around a solar- 
type st ar properly published as a "Jupiter-mass companion" 
was bv lMavor & Oueloa (11995b . T his companion to 51 Pe - 
gasi was immediately confirmed bv lMarcv & Butlen d 1995b . 
51 Pegasi b constitutes the prototype of the hot Jupiters, 
but it can also be more generally regarded as the prototype 
for all extrasolar planets discovered with the radial veloc- 
ity method. While the radial velocity method still is to be 
credited with the top score in terms of number of planets 
(and multiple planetary systems) discovered, other discov- 
ery methods gain their importance from the fact that the ob- 
servational biases involved can be significantly different. 

Direct imaging, for instan ce, is obviously biased towards 
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large semi-major ax es (e.g.|K alas et al. 2008, M arois et al. 
2008, and probably [Lagrange et al J 120091) . whereas radial 



velocity and transit measurements are biased towards the 
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Fig. 1 The "Perryman tree": Detection method s for extra-solar planets with detectable masses on a (logarithmic) mass 
scale. Adopted from Figure 1 in IPerrymanl (120001) . updated to include recent detections up to February 2010 (courtesy of 
M. A. C. Perryman). 



detection of planets on orbits with small semi-major axes 
(and large masses resp. radii). This correspondence between 
the radial velocity and transit methods is in a sense a good 
thing, since transit detections always need to be confirmed 
by radial velocity measurements in order to secure a planet 
discovery, while planet candidates from radial velocities de- 
tections usually remain candidates as long as the inclina- 
tion cannot be constrained. For ground-based surveys, the 
micro-lensing technique is most s ensitive in th e vicinity of 
the Einstein radius at 2-3 AU ([Bennett et al.ll2009t actual 
detections exist for semi-major axes in the range of 0.6 - 
5.1 AU). 

The timing method, or more precisely, the various tim- 
ing methods, also constitute indirect methods. While ex- 
ploiting the same stellar "wobbling" effect induced by an 
unseen companion that is put to use in the radial velocity 
method, it measures the varying light travel time from the 
star(s) to the observer in the course of a "wobbling" cycle. 
Its amplitude (a direct measure of the projected semi-major 
axis of the orbit followed by the central object) increases 
with large companion masses, but also with large separa- 
tions. Despite this increasing sensitivity towards wider or- 
bits, the fact that the observational time base required for 
a detection increases for longer orbital periods also must 
be factored into the overall detection probability. Exactly 
as in the case of detections from radial velocity variations, 



candidates discovered with any of the timing methods suf- 
fer from a systematic uncertainty in the mass determination 
whenever the orbital inclination remains unknown. 

A more detailed overview of the available planet detec- 
tion methods with an emphasis on their respective sensitiv- 
ity to masses has been given by IPerrymanl (120001) : an up- 
dated version with planet detection counts up to early 2010 
is given in Fig.Q] 



2.2 Planetary systems around evolved stars 

The majority of extrasolar planets known today were found 
around solar-like stars. On the one hand, this has practical 
reasons - both the number of lines in the optical and their 
sharpness decreases towards earlier spectral types, while 
later spectral types tend to show augmented spectral vari- 
ability due to activity, making small periodic radial veloc- 
ity signals harder to detect. On the other hand, any ded- 
icated searches for solar system analogues and Earth-like 
planets will obviously primarily target solar-like stars. In 
the quest to understand how our own solar system including 
the Earth has formed, it seems plausible to assume that both 
the knowledge of the frequency of similar systems as well 
as an overview of just how differently planet formation has 
proceeded elsewhere are important ingredients. 
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For stars evolving off the main sequence, the radial ve- 
locity method remains applicable in the red giant regime. It 
has turned up a total of 27 detections around G and K giants 
that add to the diversity of systems kn own. Amon g the ini - 
tial discove ries wer e t hose by Frink et al. (|2002), 



Hatzes et all d2005l) . lHatzes et all d2006l) . and 



Dollinger et al.l ( 12007b . For all stars evolved beyond the first 



red giant branch, however, the method that has most suc- 
cessfully been applied to detect extrasolar planets so far is 
the timing method. 

2.3 Timing methods 

The timing method actually comes in a variety of flavours. 
In all cases, a mechanism intrinsic to the central object pro- 
vides a stable clock the time signal of which reaches the ob- 
server with a delay or in advance to the mean arrival times 
when a further body causes a cyclic displacement (wobble) 
of the "clock". As a side note, the associated change in pe- 
riod of the clock due to the Doppler effect is typically much 
smaller than the light travel time delays. 

The most prominent example for a central object's clock 
are the pulses from rapidly rotating neutron stars. Pulsar 
pl anets were discover e d in th is way to orbit PS R 1257+12 



by W olszczan & Fra il ( 1992), and confirmed bv lWolszczan 



( 199 4J). A further d etection was reported for P SR B 1620-26 
by iBackerl (fl993l) . iThorsett et all (1 1993b . and iBacker eTal 
(1993). This system in the globular cluster M4 is known for 
its assortment of compo nents, with a white d warf in a tight 
orbit ar ound the pulsar (IThorsett et al.lll999t imaged using 
HST bv lRicher et aU2003 ), so the planet may be referred to 
as a circumbinary planet^ 

In a further approach, the timing method specifically tar- 
gets circumbinary planets by design. The clock in this case 
are frequent, sharp eclipses in a close central binary sys- 
tem. Eclipse timing has uncovered sub-stellar companions 
to subdwarf B stars (sdB, see section [4]i in HW Vir-like bi- 
naries, in (pre-)cataclysmic variables and possibly also in 
W UMa systems. Two planetary companions have been ac- 
cepted as confirmed around the prototype system HW Vir 
(ILee et all 120091) . wh ile the tertiary com ponent to the bi- 
nary HS 0705+6700 dOian et all l2009dh probably lies in 
the mass range for brown dwar fs. Further HW Vir systems 
with preliminary detections by lOian et all (120 1 Och include 
HS 2231+2441 and NSVS 14256825. For the (pre-)cata- 
clysmic variable central systems, t ertiary detections have 
been reported for the polar DP Leo ( Oian et all 2010al) . for 
the DA+dme binary QS Vir (Oianet al., 2010b), and per- 
haps for the nova NN Ser (lOian et al.ll2009ai status of the 
planet unconfirmed) |3 

1 As a side note, a different example for a hierarchical system is the 
planet in the binary 7 Cephei AB, where the planet orbits the primary K 
su bgiant component (earl y speculations from radial velocity measurements 
by ICampbell et al]|l988i were first confirmed by Hatzes et al. 20031 and 
refined through a direct detection of the secondary M dwarf component by 
iNeuhauser et al.l2007l) . 

2 Provisional reports on more detections are also out on the W UMa- 
type eclipsing binaries NY Lyr and DD Mon iOian et a l.. 2009b c). 



A variant of this method is to look for timing residuals 
in known planetary transits in order to uncover additional 
planets in the system. 

The third possibility is to resort to stellar oscillations as 
a clock. This lead to the first discover y of a planet aroun d 
a subdwarf B star, V391 Pegasi, by ISilvotti et all d2007l) . 
The existence of this system suggested the possibility that a 
planet in an orbit similar to that of the Earth may have sur- 
vived the red giant expansion of its presumably single host 
star, and has triggered follow-up searches for comparable 
systems. I will come back to this application in section [331 
and in more detail in section [5] 

3 Stellar oscillation - extrasolar planet links 



An obvious benefit of high-precision photometric time se- 
ries of planet-hosting stars (such as those obtained in ambi- 
tious transit searches) is that the host star can additionally 
be characterized in detail with asteroseismic methods if it 
oscillates. 

3.1 Asteroseismology 

For small perturbations to a spherical equilibrium solution, 
an infinite series of non-radial modes, with the radial modes 
included as a special case, can occur (which can be de- 
scribed by nodal planes in two angular directions, and spher- 
ical nodal surfaces in the radial direction), leading to multi- 
periodic frequency spectra. The equilibrium is restored by 
the actions of pressure and buoyancy, with one of the two 
usually dominating in a particular region of the star. The 
number of actually measurable frequencies depends on 
whether the corresponding modes are excited (i.e., if a driv- 
ing mechanism is converting radiative energy or, temporar- 
ily, convective movement into pulsational kinetic energy), 
and on whether they are observable as photometric or ra- 
dial velocity variations (i.e., if the mode geometry yields a 
detectable net effect in the integrated observables). 

Given appropriate modelling capabilities, the density 
structure of a star can be inferred from observations of a 
sufficiently large number of excited eigenmodes that probe 
the interior conditions differentially. While the exact anal- 
ysis approach can vary depending on the class of variables 
considered, important fundamental parameters that can be 
derived from this exercise are the stellar mass, radius and 
(depth-dependent) chemical composition. 

In solar-like pulsators, oscillatory eigenmodes are ex- 
cited stochastically by convection. The known pulsations 
are p modes (acoustic modes, with displaced material re- 
stored by the action of pressure), although g modes (grav- 
ity modes, action of buoyancy) are also thought to exist in 
the deep interior of the sun. Considering that derivatives of 
the local equilibrium gravity and radius can be neglected in 
this case, approximations for high radial order can be found. 
For acoustic modes described within the asymptotic theory, 
one finds a regular frequency spacing for modes of the same 
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low angular degree corresponding to subsequent radial or- 
ders. The large frequency separations, along with modifica- 
tions introduced by considering the effect of different de- 
grees leading to additional small frequency separations, to- 
gether allow the definition of valuable diagnostic tools. This 
includes the famous echelle diagrams, as well as the 'aster- 
oseismic HR diagram' which allows to relate the two ob- 
servables large and small frequency separation directly to 
the masses and ages of solar-like stars. 

3.2 Connection with solar-like oscillations 

The precise stellar parameters available through asteroseis- 
mic investigations contribute to answering a number of key 
questions in extrasolar planet research. One of the early 
noteworthy examples was the idea to investigate the scenar- 
ios for the origin of the enhanced metal content of p, Arae 
(HD 160691), host to a system of at least four planets, with 
asteroseismic methods. From grou n d-based radial vel ocity 
observations dBouchv et all 120051) . iBazot et al.l (12005b at- 
tempted to decide whether the overmetallicity was limited 
to the outer layers (accretion scenario) or was present 
throughout the star (resulting from enhanced metallicity in 
the original proto-stellar cloud). The latter scenario with en- 
hanced planet formation rates in intrinsically metal-richer 
star-forming clouds has now been generally accepted, a re- 
s ult incorporated into the more recent analysis by 



Soriano & Vauclairl (120091) that has in addition uncovered a 
high helium abundance. It will remain to be seen if contribu- 
tions from asteroseismology will also be able to help solve 
the possibly related mys tery of enhanced lith ium depletion 
in planet-hosting stars (llsraelian et alJ.120091) . Fundamental 
links between the two programmes of the space mission 
Corot, the exoplanet search and the asteroseismology pro- 
gramme, have been pointed out by 
Vauclair et al.l (120061) (see also lSoriano et alj2007tlVauclair 
2008). In Corot these two programs are conducted with sep- 
arate pairs of detectors operated with different instrumental 
setups: an on-focus setup with dispersion through a bi-prism 
sampled every 5 12 s in the exoplanet field, and a highly out- 
of-focus setup sampled every 1 s in the seismology field. 

While a smaller technical dichotomy continues to ex- 
ist with the long (30min) and short (1 min) cadence read- 
out modes for the individual apertures on the Kepler satel- 
lite's detectors, the boundaries start to dissolve here since 
the sampling is at best loosely associated with the classifi- 
cation of a target as belonging to the planet-hunting core 
program or a program such as the asteroseismic investi- 
gation, and can simply be re-assigned. Kepler has demon- 
strated its capabilities early on in the mission schedule with 
o bservations of the kno wn transiting exoplanet host HAT-P- 
7 (IBorucki et all 120091) . HAT-P-7 has in the meantime been 
further characterized through an analysis of its simultane- 
ously discovered solar-like oscillations 



dChristensen-Dalsgaard et all 1201 Oft 



tine analysis of solar-like oscillations in newly-discovered 
planet-hosting stars, as well as the analysis of a large variety 
of asteroseismology targets of interest to the pulsation sci- 
ence community, have been institutionalized in the Kepler 
Asteroseismic Investigation (KAI) and the Kepler Astero- 
seismic Science Consorti um (KASC, see 

Christensen-Dalsgaard et alj|2 008, and in particular the first 

1201 (J 



results 



Gilliland et al. 



_ _ . and 

Christensen-Dalsgaard et alJl2010» . The importance of de- 
termining the stellar radius as accurately as possible stems 
from the circumstance that a transit light curve yields the 
radius ratio between star and planet. In order to derive the 
transiting planet's absolute radius the value of the stellar ra- 
dius must be known. Together with the planet mass deter- 
mined from the confirmation radial velocity curve (again, 
the stellar mass must be known), the planet's mean density 
is found. In the proposed Plato mission 



(e.g. ICatala & ESA Plato Science Study Team! [2009), 



Kepler has now also de livered its first five g enuine ex- 



trasolar planet discoveries (IBorucki et all 12010). The rou 



Lgei 
li- 



the high-precision determination of planet host star radii and 
other fundamental stellar parameters from asteroseismology 
is an integral part of the planet hunting and characterization 
concept. 

3.3 Connection with coherent oscillations 

The above considerations could in principle with the same 
rationale be extended to planet-hosting stars that exhibit 
pulsations other than solar-like. Well beyond this, the strik- 
ing capabilities in particular of Kepler evidently open many 
possibilities for genuine asteroseismological applications 
that target questions in many areas of stellar astrophysics. 
This includes the classical pulsators which, instead of being 
stochastically excited as the solar-like stars, exhibit unstable 
modes driven by the k mecha nism, with t opical applications 
extensively described by Gilli land et al.l J2010l) . 

The standard asteroseismic exercise derives the instan- 
taneous structure of a star, often relying on model structures 
from full evolutionary calculations in the process. As an ex- 
tension, period changes in the oscillatory eigenfrequencies 
due to evolutionary effects can be considered as an addi- 
tional constraint to find the best solution in parameter space. 
Given a series of evolutionary models already subjected to a 
stability analysis, the rate of change P of modes in a specific 
model can be determined without too much trouble. Due to 
the typically very long evolutionary time scales involved, 
measuring the secular evolution in real stars is observation- 
ally expensive and somewhat more complex. 

An example of a class of objects where pulsators can be 
found that prove to be coherent and stable on time scales of 
many years are the ZZ Ceti variables on the DA white dwarf 
cooling track. The pulsations in these objects are due to the 
recombination of hydrogen in a narrow temperature range, 
leading, via the associated increase of the opacity in the 
outer layers, to the manifestation of low-degree low-order g 
modes. In contrast to the millisecond pulses in spun-up pul- 
sars, and also in comparison to the pulse duration in normal 
"slow" pulsars, the pulsation periods in white dwarfs are 
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much longer: of the order of a few minutes. The precision 
that can be reached in determining the clock rate is hence 
accordingly lower, and long-term changes are more readily 
analyzed using O— C techniques, instead of determining P 
directly as the derivative of a series of quasi-instantaneous 
P measurements over time. 

Measurements of P and its interpretation in the con- 
text of cooling times exist for a small number of suitable 
pulsating white dwarf s: PG 1 159-035 dCosta et alj|l999 : 



Costa & Kepler! 120081). Gl 17-B15A J Kepler et all 



1991 



Kepl er et alj 12000: Kep ler et al.l 12005). and a larger sam 



pie of a total of 15 objects investigated by Mullally et al 
(120081) . It was quickly recognized (e.g. by IProvencall 1 1 9971) 



that the influence of a possible unseen companion can be 
measured as a side effect from the same data. Applying the 
ideas of the timing method ( section [23b has given the long- 
term photometric monitoring of pulsations in white dwarfs 
and related objects a new spin as a means to search for plan- 
ets aroun d evolved stars. With respect to previous work, the 
efforts bv lMullallv et all (2008) show a shift of focus to that 
effect. Their proposed planet candidate around GD 66 has 
however remained unconfirmed. 

The same measurements are possible for a different 
group of compact oscillators, the subdwarf B stars (see be- 
low), a quantitatively uncommon feeder channel for white 
dwarfs. The rapid pulsations in subdwarf B stars, of the 
order of minutes just as in the ZZ Ceti white dwarfs, are 
generated via a n mechanism, providing potentially suitable 
conditions for reasonable long-term coherence. Yet, results 
have only been published for one pulsating subdwarf B star 
so far, V391 Pegasi (aka HS 2201+2610). As stated in sec- 
tion l2.3l these measurements, simultaneously to P, revealed 
the presence of a companion with a planetary mass. The 
possible reasons for this initially unexpected, instant suc- 
cess, sometimes jokingly referred to as "100% discovery 
rate", are worth a more in-depth investigation. 

4 Hot subdwarf stars 

4.1 Evolution 

Subdwarf B stars (sdBs) are subluminous hot stars that are 
found in an effective temperature range from 20 000 K to 
40 000 K at surface gravities between about 5.0 and 6.2 in 
log (<?/cms~ 2 ), and that can in many cases be identified 
with evolved stellar models on the extreme horizontal 
branch (EHB). Their masses are expected to peak around 
the value for the He core flash at 0.46 M Q . As is true for 
all horizontal branch stars, extreme horizontal branch stars 
have a He burning core but, due to previous significant mass 
loss, no H-shell burning in their thin hydrogen shells. The 
thinness of the shell leads to their blue appearance, so that 
the flux from hot subdwarfs (including the hotter sdOs) con- 
tributes significantly to the UV excess observed in galaxy 
bulges and elliptical galaxies. A great recent review of the 
observational properties of hot subdwarfs including bina- 
rity, kinematics, as well as current modelling capabilities 



to describe their atmos pheres, interio rs, and evolutionary 
history, can be found in IHeberl d2009l) . The most puzzling 
question about hot subdwarfs remains what the precise evo- 
lutionary status of these core-helium burning or even more 
evolved objects really is. Here I only focus on how to pos- 
sibly produce the sdB type. The basic problem in the stan- 
dard single-star evolutionary scenario is to shed the hydro- 
gen envelope almost entirely just before or at the moment 
of the He flash. This would require ad hoc strong mass 
loss through stellar winds for a certain fraction of stars upon 
reaching the tip of the first giant branch. Explaining the for- 
mation of subdwarf B stars in the context of binary evo- 
lution is t herefore now g enerally favoured over single-star 



scenarios (Ha n et aUl20 03, 2002) 



Following the overview bv lPodsiadlowski et alj ([2008), 
three genuine binary formation scenarios, and additionally 
the merger scenario, can be distinguished. All three of the 
binary scenarios involve Roche-Lobe overflow 
(RLOF), each at different stages and under different con- 
ditions. 

In the "stable RLOF + CE" channel, the system initially 
goes through a first mass-transfer phase with stable RLOF 
that turns the evolving component into a He white dwarf. 
When the second component, the future subdwarf B star, 
then evolves to become a red giant, the second mass-transfer 
phase can happen dynamically. This unstable RLOF, where 
the matter transferred cannot all be accreted by the He white 
dwarf, leads to the formation of a common envelope (CE). 
After spiral-in and envelope ejection, the resulting system 
consists of the He white dwarf and the sdB - the core of the 
giant with its envelope removed - in a short-period binary 
(< 10 days). 

In the "CE only" channel, unstable RLOF occurs when 
the future subdwarf B star starts transferring matter to a 
lower-mass main-sequence companion near the tip of the 
first giant branch. The ensuing common-envelope phase 
again leads to a closer final configuration of the resulting 
system, which will consist of a low-mass main sequence 
star and the sdB in a short-period binary (< 10 days). 

Besides these two common-envelope channels (involv- 
ing unstable RLOF), the "stable RLOF" channel can also 
produce sdBs. Stable Roche-Lobe overflow can occur when 
the future subdwarf B star starts transferring matter to a 
main sequence companion at mass ratio below ~ 1.2. As 
before, this happens near the tip of the red giant branch, 
but this time the system more likely widens due to the mass 
transfer. The resulting system consists of a main sequence 
or subgiant star and the sdB in a wide binary (;> 10 days). 

The problem with all of the above binary scenarios is 
that (apparently?) single sdBs also exist. In addition to the 
single star scenario involving variable mass loss as men- 
tioned above, a further possibility to produce these are 
mergers. In the "merger" scenario, two He white dwarfs 
in a close system, produced while undergoing one or two 
common envelope phases, spiral towards each other due to 
angular momentum loss via gravitational radiation until the 
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less massive one gets disrupted and its matter accreted onto 
the more massive component. At a critical mass, the accre- 
tor can ignite helium fusion and the merger product would 
hence indeed turn into a single sdB. 

In the confrontation with observations, both the types 
of companion found in close systems (white dwarf or low- 
mass main-sequence stars) as well as the close binary fre- 
quency are roughly as expected when the above merger 
channel is considered. Yet observationally, the mass spec- 
trum for the companions is broader than expected from the 
standard formation scenarios for sdB+WD, sdB+dM and 
single sdBs. On the low-mass end, it can be argued to com- 
prise the V391 Pegasi planet and, possibly, further 
secondary sub-stellar companions. Four companions with 
unusually h i gh masses have als o been reported; according to 
Geier et al.l (12009a. l2010l 120081) . the massive compact com- 
panions found from radial velocity variations must at least 
be heavy white dwarfs or in two cases even neutron stars or 
black holes. 

While the problem of sdB formation has not been fully 
solved, the subsequent evolution of a subdwarf B star to- 
wards the white dwarf cooling sequence is more straight- 
forward. One interesting characteristic is that this evolution, 
due to the thinness of the outer layers, bypas ses the asymp 



A number of objects are known that show both types of 
mode simultaneously and are referre d to as hybrid pu l sators , 



totic g iant branch, as was noted early on by Dorman et al 



(119931) . Another aspect is that the sdB+dM close systems 
constitute potential progenitors for cataclysmic variables. 
Furthermore, given the variety of configurations in which 
sdBs are found as one component, the diversity of possible 
progeny systems is not restricted to pre-CVs. 

All of the above, while far from understood in the de- 
tails, implies that the sdB's stellar core as it looks like after 
the first giant branch - or after a merger - is almost laid 
bare. It is hence very well accessible to asteroseismological 
methods. 

4.2 Pulsating subdwarf B stars 

Only a small fraction of the sdBs show pulsational varia- 
tions, with non-pulsators also populating the region in the 
HRD where the pulsators are found. There are p- (pressure-) 
mode and g- (gravity-)mode types of pulsation. 

The rapidly pulsating subdw arf B stars (sdBVr) were 
discovered observationally by Ki lkenny et al.l (119971) . The 



short periods of these p-mode pulsators are of the order of 
minutes and have amplitudes of a few tens mmag. These 
pulsations were inde pendently predicted to exist by 



Charpinet et al.l (119971) . who in this initial and subsequent 
research papers explain the driving to result from a k mech- 
anism due to a Z-opacity bump accumulated by radiative 
levitation. 

The longer periods in the slowly pulsating subdwarf B 
stars (sdB Vs) range from 30 to 80 mins at even lower ampli- 
tudes of a fe w mmag. Thi s group of g-mode pulsators was 
discovered bv lGreen et al . (2003) and has been explained to 
pulsate due to t he sam e k mechanism as the sdBV r type by 



or sdBVrs in the nomenclature of iKilkenny et al.l (120101) 



that is used here. These hybrids lie in the overlapping re- 
gion between the hotter short-period and the cooler long- 
period pulsators. The prototype and a prominent example 
is HS 0702+6043. It s rapid pulsations were discovered by 



Dreizler et al. (2002), and subsequently the slower ones by 
Schu h et al.l ( 2006, 2005). Further examples can be found in 
Schuhl d2008l) and references therein. It should however be 
noted here that the planet host star V391 Pegasi 
(HS 22 01+2610) is also among the hybrids. It ha s five p 
modes dOstensen et al.ll20 01b; Silvotti et al. 2002a) and at 
least one g mode (lLutz et alll2009ai l2008a). 

Exploiting the sdBVr asteroseismologically has the po- 
tential to test diffusion processes and has lead to a first mass 
distri bution for subluminous B stars (ICharpinet et al.ll2008 ; 
Fontaine et al. 2008; Ostensen 2009 and references therein.) 
On the other hand, the extent of the instability region for the 
sdBVs, but in particular the existence of the sdBVrs, have 
challenged details of the input physics for the models. The 
actual composition of the 'Z' in the Z bump (iron, nickel) 
as well as the role of opacities hav e been discussed in this 
context dJefferv & Said. 120061 12007b . While, quite similarly 
to the findings for the sdB variables, using updated Opacity 
Project (OP) instead of OPAL opacities also improves the 
situation in (3 Cep and [SPB] variables, the opposite is the 
case in the sun with helioseismology. The most appropriate 
opacities are therefore still under discussion. 

Further open questions in this field are the co-existence 
of pulsators and non-pulsators, and also the origin of the 
amplitude varia bility observed in a number of sdBVs 
dKilkennvil2O10h . 



4.3 Sub-stellar companions of subdwarf B stars 

When a companion with planetary mass was found around 
the hybrid pulsating subdwarf B star V391 Pegasi with the 
timing method, this indicated that a previously undiscovered 
population of sub-stellar co mpanions to apparently sing le 
subdwarf B stars might exist (iSilvotti et all l2002al 120071) . 

As for low-mass tertiary bodies aro und HW Vir -like sys- 
tems detected by eclipse timing, iLee et al.l d2009l) have re- 
ported two p lanetary companio ns around the system 
HW Vir itself. lOian et all (|2009d) put the mass of the ter- 
tiary body in the HS 0705+6700 system in the range for 



brown dwarfs. Still unconfirmed detections by Qian et al. 



d2010d) exist for HS 2231+2441 and NSVS 14256825. 
A different type of d etection has been put forward by 



Fontaine etal .1 (12003) 



Geier et al.l d2009a[ l2010» . who report measurements of ra- 
dial velocity variations in HD 149382 indicative of a plan- 
etary companion, variations which if real are indicative of 
a close sub-stellar companion in a 2.4 d orbit. Based on 
the above and this finding (that yet has to be confirmed), 
these authors also argue for a decisive influence of sub- 
stellar companions on the late stages of stellar evolution. 
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Fig. 2 Photometric data obtained with the 1.2 m MONET/North telescope during a period of 2 weeks in May 2005, 
showing the pulsations in the light curve of V391 Pegasi. Starting from a mean intensity of unity on the first night, subse- 
quent observations are shifted downward by a fixed offset each night. The actual data points (black dots) are overlaid with 
model for the pulsations in red. One such run spanning several nights is required to derive one O— C point by comparing 
the current phasing to that of a mean model. 



The EXOTIME program aims to increase the empirical data 
available on which to base such discussions. 



5 The EXOTIME program 

5.1 The planet-hosting pulsating sdB V391 Pegasi 

V391 Pegasi (HS 2201+2610) was first discovered to be a 
rapidly pulsating subdwarfB star by |0stensen et al.l ( 12001a) . 
Additional slow p u lsations were subsequently r eported by 
Lutz et alJ d2009ah . ISilvotti et all (J2002al 120071) were able 
to derive P values for the two strongest pulsation modes, 
and found an additional pattern in the observed-calculated 
(O— C) diagrams that revealed the presence of a giant planet 
in a 3.2 year orbit. The fact that this cyclic variation has 
been measured independently from two frequencies consid- 
erably streng t hens th e credibility of this discovery. Actually, 
Silvotti et al.l ((2007) detected parabolic and sinusoidal vari- 
ations in the O— C diagram constructed for the two main 
pulsation frequencies at 349.5 s and 354.2 s over the ob- 
serving period of seven years. The sinusoidal component 
with its 3.2 year periodicity is attributed to the presence of 
the very low-mass companion V391 Pegasi b at msini = 

3.2 ± 0.7Mj u _. The scenarios proposed for the origin of 
this planet are discussed in section[i 



5.2 Further characterization of V391 Pegasi system 

The determination of the true mass of the "asteroseismic 
planet" V391 Pegasi b requires to find a constraint on the or- 
bital inclination. Besides the orbital inclination, the orbit ec- 
centricity has also not been well determined so far. The con- 
straints on possible further planets are weak and currently 
allow for a second planet in the system massive enough to 
be detected. Continued photometric monitoring (on-going, 
see for example Fig. 12 Silvotti et al. in prep.) will be able 
to: check if the O— C evolves as predicted from the orbital 
solution, investigate the eccentricity, investigate a possible 
multiplicity, attempt another independent re-detection of the 
planet from P3, and search for rotational splitting in the pul- 
sations (see below). 

As a first step to derive the mass of the known compan- 
ion object with follow-up observations. ISchuh et alJ (120091) 
have attempted to determine the orbital inclination from 
spectroscopy. The approach suggested to achieve this is to 
use the stellar inclination as a primer for the orbital orien- 
tation. "Stellar inclination" can refer to the rotational or the 
pulsational axis, which as a further necessary simplification 
are assumed to be aligned, and can in turn then be derived 
by combining measurements of v to ^ and v ro i sin i. 
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Fig. 3 Detail of the observed Ha core in a pulsation-integrated Keck spectrum of V391 Pegasi. The data (shown as 
a black line) are plotted centered on 6562. 798A. While measurements on the Ha core may be subject to systematic 
effects, this line is the only one available for such investigations in the spectral range covered by the observation. The 
flux was normalized; it does not reach unity at the edges of the plot as these regions are still well within the line wings. 
Overlaid are 1) a model spectrum in red calculated from a NLTE pure H atmosphere model at T e ff = 30 000 K and 
log (g/cms~ 2 ) = 5.50, and 2) a model spectrum in blue from the same atmosphere model as before but rotationally 
broadened with . 2 A corresponding to 9 km s ~~ . 



The value for v ro i is in principle accessible through ro- 
tational splitting in the photometric frequency spectrum 
(which has however not been found for V391 Pegasi yet), 
while the projected rotational velocity v ro i sin i can be mea- 
sured from the rotational broadening of spectral lines. This 
rotational broadening must be deconvolved from the addi- 
tional pulsational broadening caused by the surface radial 
velocity variation in high S/N phase averaged spectra. 

Both phase averaged and phase resolved high resolution 
echelle spectra were obtained in May and September 2007 
with the Hobby-Eberly Telescope (HET), and one phase av- 
eraged spectrum in May 2008 with the Keck 1 telescope, 
in order to put limits on the pulsational radial velocities. 
Echelle spectra of V391 Pegasi were taken during May and 
September 2007 with the HRS (R = 15 000) of the HET at 
the McDonald Observatory, and with HIRES (R = 31 000) 
at the Keck 1 telescope atop Mauna Kea in May 2008. 

Using standard data reduction procedures, the individ- 
ual echelle orders were merged and the final spectra care- 
fully normalized and finally summed. This results in a set 
of individual spectra (S/N « 3), in particular two Septem- 
ber 2007 high time resolution series, and summe d spectra 
for May and September 2007 (tKruspe et alll2008l) . 

In an attempt to "clean" the relevant rotational broaden- 
ing from pulsational effects, the spectra in September ob- 
tained in time resolved mode were combined to a series of 
ten phase resolved averaged spectra ( S/N ss 9) for the m ain 
pulsation period of 349.5 s (similar to lTillich et al.ll2007l) . 

The cross-correlation of this series of averaged spec- 
tra with a pure hydrogen NLTE model spectrum at T e ff = 
30 000 K and log (g/cms~ 2 ) = 5.5 as a template yields 



pulsational radial velocity measurements for the different 
Balmer lines. The maximum amplitude of a sinusoidal curve 
(fixed at the expected period) that could be accommodated 
in comparison to the weighted means of the Balmer lines re- 
veals that any pulsational radial velocity amplitude is smal- 
ler than the accuracy of our measurements and confirms the 
upper limit of 16 km s _1 given bv lKruspe et alj ([2008). 

The resolution of the model template matches the spec- 
tral resolution of the (pulsation-averaged) Keck spectrum. 
A comparison of the Ha NLTE line core shape (see Fig. [3]) 
yields a more stringent upper limit for the combined broad- 
ening effect of pulsation and rotation of at most 9kms~ 1 , 
meaning better spectral resolution and signal to noise data 
will be necessary to measure f DU ] S and v TO i sin i. 

5.3 Exoplanet Search with the Timing Method 

Following the serendipitous discovery of V391 Pegasi b, 
the EXOTIMBj monitoring program was set up to search 
for similar systems. EXOTIME monitors the pulsations of a 
number of selected rapidly pulsating subdwarf B stars on 
time-scales of several years with the immediate observa- 
tional goals of a) determining P of the pulsational periods 
P and b) searching for signatures of sub-stellar companions 
in O— C residuals due to periodic light travel time varia- 
tions, which would be tracking the central star's companion- 
induced wobble around the centre of mass. 

The long-term data sets should therefore on the one hand 
be useful to provide extra constraints for classical astero- 



http: //www. na. astro . it/~silvotti/exotime/ 
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Table 1 Overview of the EXOTIME targets, see section [5] for details. 



target 



coordinates (equinox 2000.) mg status 



HS 2201+2610 22:04:12.0 +26:25:07 

HS 0702+6043 07:07:09.8 +60:38:50 

HS 0444+0458 04:47:18.6 +05:03:35 

EC09582-1137 10:00:41.8 -11:51:35 

PG 1325+101 13:27:48.6 +09:54:52 



14.3 collecting data 

14.7 collecting data 
15.2 collecting data 
15.0 collecting data 

13.8 collecting data 



aka V391 Pegasi, planet candidate, sini unknown 
see Schuhetal. (2010) 
see Schuhetal. (2010) 



seismological exercises from the P (comparison with "lo- 
cal" evolutionary models), and on the other hand allow to 
investigate the preceding evolution of apparently single sdB 
targets in terms of possible "binary" evolution by extending 
the otherwise unsuccessful search for companions to poten- 
tially very low masses. 

As noted before, timing pulsations to search for com- 
panions samples a different range of orbital parameters, in- 
accessible through orbital photometric effects or the radial 
velocity method: the latter favours massive close-in com- 
panions, whereas the timing method becomes increasingly 
more sensitive towards wider separations. A further advan- 
tage of timing versus radial velocities is that the former, al- 
though observationally expensive, is easier to measure than 
the latter. In fact, it is very hard to achieve the required ac- 
curacy in radial velocity measurements from the few and 
broad lines in hot subdwarf stars (a notable exception is 
the publication of small RV variations in the bright sdB 
HD 149382 as reported bv lGeier et al.ll2009ah 

The targets selected for monitoring in the EXOTIME 
program are listed in Table [T] The target selection criteria 
ap plied to compile th is list over time have been described 
bv lSchuhet ail d2010l) . 

V391 Pegasi (HS 2201+2610) appears as the first en- 
try, as the monitoring of this system is on-going (see 15. It . 
The rapid pulsations in HS 2201+2610 were discovered by 
Ostensen et al. (2001b), and additional slow pulsations by 
Lutz et al.l (2009a). An asteroseismolog y analysis of the star 
is included in Tsilvotti et al] d2002all2007h . 

Rapid oscillations were di scovered in the secon d target 
on the list, HS 0702+6043, by breizler eTall d2002h. and si 



multaneous slow oscillations were reported by lSchuh et al. 
J2006) . The on-going EXOTIME observations for 
HS 0702+6043 have also previously been summarized by 
Lutz et al.l (1200 8b, 2009b) and include a significant contri- 
bution of data by lFrancoeur et al.l (I2010I) . 

The t arget HS 0444+0458 was first discovered to pul- 

I " 1 1 — L 

sate by Ostensen et al. (2001a), and has been further char- 

acterized bv lReed et al.1 120071) . EXOTIME has followed it 

reg ularly since 200 8. 

iKilkenny et al.l (I2006J) discovered rapid pulsations in 
EC 09582- 1137 which iRandall et all (12009b subjected to 
an asteroseismology analysis. EC 09582— 1137 is included 
in EXOTIME as a southern hemisphere target. 

The discovery of pulsations in PG 1325+101 by 



dCharpinet et all 120061 ISilvotti et all |2006ij). This target is 
also extensively being observed within EXOTIME. 

Schuh et al.l OOlOh present a portion of the observations 



currently available, describe the treatment of the data and 
display the first, still relatively short, O— C diagrams for 
the EXOTIME targets HS 0444+0458 and HS 0702+6034. 
Not surprisingly, these illustrate the need for further obser- 
vations and i nterpretation. The analysis reso rts to too ls pro- 



vided by [M ontgome ry & O'Donoghud (119991) and 



Lenz&Bregerld2005l) 



Silvotti et al.l d2002bl) was also taken advantage of by char- 
acterizing the star asteroseismologically 



6 Sub-stellar companions of evolved stars 

6.1 Post-main-sequence evolution 

For single stars with initial masses below ~ 2.2 M Q on the 
main sequence, the canonical steps in the evolution are: hy- 
drogen shell burning following core hydrogen exhaustion, 
core mass increase (due to the shell burning) and core con- 
traction with increasing electron degeneracy simultaneous 
to radius swelling on the first giant branch (RGB), ended 
by onset of a core helium flash, which brings the star onto 
the horizontal branch. Their masses now lie in the range 
roughly 0.5 M© < M < 1.0 M Q , with a canonical core 
mass of 0.46 M Q and variable envelope masses with or 
without hydrogen shell burning that determine the location 
on the horizontal branch. The hottest, bluest objects are the 
least massive and have thin envelopes that put the most ex- 
treme of these objects with inert hydrogen shells on the 
EHB. This sequence is terminated at the point where the 
helium main sequence coincides with the horizontal branch 
helium core mass. 

After core helium exhaustion the complex asymptotic 
giant branch evolution with H and He shell burning, ther- 
mal pulses, convection and dredge-up processes, and ex- 
treme mass loss unfolds. The subsequent post-AGB phase 
sees the creation of a planetary nebula from the envelope 
material lost and the emergence of the stellar remnant on its 
way to the white dwarf cooling track. 

As the subdwarf B stars are in between the RGB and 
AGB giant expansion phases, or more precisely, will not un- 
dergo AGB evolution, their observation allows to separate 
the effects of the first giant branch on planets from those of 
the asymptotic giant branch. 
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Fig. 4 The subdwarf B (sdB) star V391 Pegasi oscillates in short-period p modes and long-period g modes, making it one 
of the known hybrid pulsators among sdBs. As a by-product of the effort to measure secular period changes in the p modes 
due to evolutionary effects on a time scale of almost a decade, the O— C diagram has revealed an additional sinusoidal 
component attributed to a periodic shift in the light travel time caused by a planetary-mass companion around the sdB star 
in a 3.2 year orbit. In the above artistic impression, the V391 Pegasi system is shown at an earlier evolutionary stage in 
one of the proposed scenarios where, roughly 10 8 years ago, the star, at maximum red giant expansion, almost engulfed 
the planet. © Image courtesy of HELAS, the European Helio- and Asteroseismology Network, funded by the European 
Union under Framework Programme 6; Mark Garlick, artist. 
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6.2 The fate of planets around evolving stars 

The starting point of discussions on planets around evolved 
stars is often that relic planets around white dwarfs are ex- 
pected simply due to the considerable number of planets 
found around main sequence stars. The question then is at 
what initial masses and orbital separations the planets ac- 
tually can survive both giant expansion phases of their host 
star with out being engulfed, di srupted, evaporated or 
ejected. iDebes & Sigurdssonl d2002l) consider the effects of 
mass loss, planet-planet interactions, and orbit stability and 
conclude that while inner planets will perish, far-out small 
bodies and distant planets have the potential to create a new 
dust disk. Such as disk could then pollute the central white 
dwarf's atmosphere to create a DAZ, and be observable via 
its IR excess. - While the search for planets around single 
white dwarfs has not been successful yet, the existence of 
dusk disks has indeed been established obser vationa l ly. 



In addition to orbital effects, IVillaver & Livid (2007) 
consider the thermal conditions during the planetary nebula 
phase and establish new regions in the orbit and mass pa- 
rameter space where planets survive all the way to the white 
dwarf stage of their host star s. Focusing more on the or- 
bital evolution, IHaghighipourl (12008b investigates the issue 
of survival for planets in binary systems. A consequence in 
all of these studies is that finding planets at parameter com- 
binations that correspond to previously "forbidden regions" 
immediately requires that second generation planet forma- 
tion scenarios must have been at work. This can in principle 
include migration and significant accretion of previously ex- 
isting planets in newly formed disks. 

Of course, these investigations are predominantly con- 
cerned with the more violent effects during the AGB phase. 
However, once a planet has made it safely into the orbit of a 
subdwarf B host star, it can be expected to continue to exist 
without too much hassle when the host evolves to its final 
white dwarf stage. 

A first scenario that explained the case of V39 1 Pegasi b 
basically in the context of single star evolution assumes that 
the planet is an old first generation planet which survived 
a common envelope phase. The situation where the star at 
maximum red giant expansion almost engulfed the planet is 
depicted in Fig. [4] Assuming this scenario, the existence of 
V391 Pegasi b proves observationally that gas-giant planets 
can survive the first red giant expansion in orbits similar to 
that of our Earth. 



Silvotti (2008) puts forward alternative scenarios. In a 



second variant, V391 Pegasi b would have been an outer 
planet, never really under a strong direct influence of the 
expanding host star. The exceptional mass loss of the host 
star, turning it into a sdB, would instead have been triggered 
by another closer-in planet that got potentially destroyed in 
the process. 

Scenario three explains V391 Pegasi b as a young sec- 
ond generation planet formed in a disk resulting from the 
merger of two He white dwar fs. The planet w ould then be a 
helium planet (compare also lLivio et al.ll2005l) . 



In the case of the planets around HW Vir-like systems, 
the situation is somewhat easier as there is no dispute about 
the origin of the enhanced mass loss in a c ommon envelope 
ejection phase. Besides the argument by iLee et alj (120091) 
that the planets were formed in a circumbinary disk, this 
leaves room for an interesting variety of possible planet for- 
mation scenarios, including a second generation scenario 
(cited from Heber 2009): "At birth, HW Vir's binary components 
must have been much further apart than they are today. During 
this [red giant] mass-loss episode the planets as well as the stellar 
companion may have gained mass. Rauch (2000) suggested that 
the low-mass companion in AA Dor may have grown from an initial 
planet by mass accretion during the CE-ejection phase. Could that 
have happened to the HW Vir system? As a matter of pure spec- 
ulation, the star may have been born with three massive planets: 
the innermost launched a CE ejection, spiraled in, and accreted so 
much material as to turn into a low-mass star. [. ..] Another specu- 
lation concerns in situ formation; that is, could the planet have been 
formed during the CE phase?" 

6.3 The fate of evolving stars with planets 

As already implied above, it has recently been suggested 
that the presence of planets may have implications for the 
formation of sdB stars. This marks the renaissance of an 
older idea: Soker (1998) suggested that planets may con- 
stitute the second parameter influencing the irregular mor- 
phology of the (blue part of the) horizontal branch in glob- 
ular clusters and elliptical galaxies. The physical processes 
discussed as candidates for the second parameter in hori- 
zontal branch morphology are 

- age of the globular cluster, 

- deep helium mixing and radiative levitation, 

- large He abundance or fast rotation, 

- stellar density in the cluster, 

- planets enhancing mass loss on the red giant branch. 
Quite similarly to this last suggestion dSoken, Il998l) . ISoker 
( 120 101) points out that in addition to strengthening the blue 
component of the normal horizontal branch, enhanced mass 
loss due to planets could also be decisive in forming sub- 
dwarf B stars. The initial prediction - that the engulfed 
planet that helps to shed the envelope will in most cases be 
destroyed in the process - remains valid, so this hypothesis 
is challenging to test. 

However, V391 Pegasi, while most probably not capa- 
ble to have caused the enhanced mass loss of its host star it- 
self, may well be regarded as a possible tracer for (former?) 
inner planets that could indeed have been able to cause the 
amount of mass loss required (compare the second scenario 
in the previous section). The scenario as a general expla- 
nation for the formation of a significant part of the single 
sdB star population may have a number of unsolved prob- 
lems, mostly related to our poor understanding of the com- 
mon envelope ejection mechanism. But, independently of 
the theory, EXOTIME and other programs now start to built 
up an empirical data base of low-mass companion statis- 
tics that further discussions can be based on. This should 
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help to clarify the role of an - to date perhaps mostly still 
undiscovered - population of very low-mass companions to 
apparently single subdwarf B stars in the formation of these 
objects. 

Should it one day become possible to determine the 
composition of the low-mass objects found, and should 
these turn out to be planets that were most likely formed 
in the helium disk of a He white dwarf merger event, the 
origin of the host stars will also have been cleared up unam- 
biguously. 

7 Summary 

This article has attempted to highlight the added value of 
incorporating stellar pulsations in the comprehensive inves- 
tigation of various stellar and planetary systems. 

The potential of links between asteroseismology and 
exoplanet science was first illustrated for solar-like oscilla- 
tors. Accurate asteroseismic radii for (solar-like) host stars 
of transiting exoplanets translate into good radius determi- 
nations for the transit planets. 

In the context of connections with coherent pulsations, 
as found in evolved stars, the focus was directed towards the 
investigation of stellar and planetary systems at late stages 
of stellar evolution. The asteroseismology-exoplanet-con- 
nection was explored in particular for post-RGB stars. The 
exotic class of subdwarf B stars, of which many show pulsa- 
tions, was introduced as a group with several new planet dis- 
coveries. In discussing the fate of planetary systems based 
on known planets around extreme horizontal branch stars 
(e.g. V391 Pegasi), the question arises whether these are 
first or second generation planets. Data from programs such 
as EXOTIME, which make use of the timing method to find 
more such systems, can fill a gap with respect to the sensi- 
tivity to planets in wide orbits and applicability to evolved 
stars. 

With respect to the puzzle of subdwarf B formation, 
planets were proposed to (partly) be at the origin of sub- 
dwarf B stars, and also to play a role as the second param- 
eter for the horizontal branch morphology in globular clus- 
ters and elliptical galaxies. 

The aspects this article summarizes are all from the pre- 
Kepler era. It will be fascinating to see how our understand- 
ing in a great many areas evolves with the new data. 
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